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A. Introduction and Overall Svunmary 

The objective of the work performed under NASA Grant NGR 22»012-031 
has been to examine both theoretically and experimentally the nature of 
the acoustic field of a simple source In a wind tunnel under flow con- 
ditions. The motivation of the study has been to seek to establish 
aspects of the theoretical framework for Interpreting acoustic data 
taken in wind tunnels using "In wind" microphones. 

There are three distinct investigations which have been performed 
under this grant. In the first which is described in Section B, the 
wind tunnel is modeled by a rectangular duct of infinite length. A 
point source is located Inside the duct, and there is a uniform flow in 
the duct. The walls of the duct are characterized by an acoustic 
impedance. A steady-state sound field Is assumed to exist in the duct. 
Owing to mathematical difficulties of this problem, the solution for 
the acoustic pressure field has not been previously calculated. An 
exact solution to this boundary-val^e problem however has been obtained 
in this work using a numerical procedure. In this solution there is no 
limitation placed on the value of the acoustic Impedance of the walls. 
Hence, the effects of both acoustically hard and soft walls may be 
investigated. Owing to computer limitations, the maximum source fre- 
quency is approximately 250 Hz for a duct of dimensions 2 meters by 3 
meters. As such this exact solution describes a "near-field" situation 
in a wind tunnel test. 

The second investigation which is described in Section C, is also 
theoretical. The acoustic pressure, acoustic power, and Doppler shift 
in frequency of a moving point source in a moving acoustic medium are 


calculated. The subsonic source and medlua speeds are arbitrary. This 
is a far^fleld calculation which neglects all wall effects. As such it 
pertains to high frequency measurements of noise using **ln wind" micro- 
phones which are placed within the "hall" radius of the wind tunnel. 

This situation may be viewed as a simple model of jet noise measur^ents 
using "in wind" microphones. In part of this work, an existing expres- 
sion for the acoustic Intensity in a uniformly moving medium has been 
shown to be Inappropriate for this free-field calculation. 

Described in Section D is the third study which is an experimental 
investigation of pulse propagation in the 2.13 meter by 3.05 meter 
(7 feet by 10 feet) wind tunnel of the NASA Ames Research Center. Low 
frequency acoustic pulses are recorded by microphones in the presence 
of flow in the wind tunnel. These data are compared vrith the results 
of a theoretical model of pulse propagation in a long duct with acous- 
tically hard walls. Both the experimental results and the theoretical 
results are in very good agreement with each other. 

Detailed accounts of these investigations either currently appear 
or will appear at a later date in various forms. Some preliminary work 
on the material described in Section B appears in an article entitled 
"Acoustic Field cf a Point Mass Source in an Infinitely Long Rectangular 
Duct with Flow" by J. E. Cole and I. I. Sarrls published in the Proceed- 
ings of the Second Interagency Symposium on University Research in 
Transportation Noise (June, 1974). An article entitled "Theoretical 
Results Pertinent to Measurements of Jet Noise in Wind Tunnels" by 
J. E. Cole and I. I. Sarris which is based on the. results of Section C 
and their relationship with jet noise measurements in wind tunnels 
appears In the Proceedings of the Third Interagency Symposium on 
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Unlverslty Research in Transportation Noise (Kovenber 1975). A manu- 
script based on Section D and entitled "Acoustic Pulse Propagation in 
a Duct with Plow" by J. E. Cole was submitted in October 1975 to the 
Journal of Sound and Vibration for publication. A letter to the editor 
based on part of the work in Section C and entitled "Acoustic power of 
a moving source in a moving medium" by J. E. Cole and I. I. Sarrls was 
submitted in October 1975 to the Journal of the AcoustJc>.l Society of 
America for publication. A detailed description of the work of Sections 
A and B will be Included in a doctoral dissertation of I. I. Sarrls who 
has been supported as a Research Assistant by this grant. It is also 
anticipated that a manuscript on the material of Section A will be 
prepared subsequently fcr publication. Finally, various aspects of 
this work were discussed in two semi-annual progress reports dated 
February 28, 1974 and October 15, 1974. 

The NASA Technical officer for this grant is Warren F. Ahtye, 
Large-Scale Aerodynamics Branch, NASA Ames Research Center. 

B. Low-to-mid Frequency Theoretical Solution 

In this section the calculation is described of the acoustic field 
of a point source located in a long rectangular duct carrying a uniform 
flow. For this calculation the walls of the duct may be characterized 
by a constant acoustic impedance. In order to investigate the acoustic 
field of the sound source, the following assumptions are made: 

(1) The medium is homogeneous, inviscid, non-heat conducting 
and in uniform motion. 

(2) The source of sound is a fluctuating point mass. 
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(3) The duct is of rectanguler cross section and of infinite 
length. 

(A) The acoustic characteristics of the walls of the duct are 
represented by a constant impedance. The two pairs of 
opposite walls may have different acoustic Impedances. 

(5) The acoustic quantities have a harmonic time dependence. 

The first two assumptions allow the linearized continuity and 
momentum equations to be combined to provide the following governing 
equation: 
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where a is the speed of sound in the medium at rest, p “ p(x,y,z,t) 
o 

is the acoustic pressure disturbance, V is the three-dimensional 
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gradient operator, ^ ^ "convective" time derivative, 

and U is the speed of the acoustic medium in the positive x-dlrectlon. 
For a harmonic point mass source located at source 

strength q is given by 


q(x,y,z,t) - 6(x-x^)6 (y-y^)fi(z-z^) 


( 2 ) 


where q^ is the mass rate, u is the source driving frequency and 
6 is the Dirac delta function. 

In accordance with the fourth assumption, the boundary conditions 
to be satisfied by the acoustic pressure at the walls of the duct are: 
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where b and d are the height and width of the duct and B » B 

y * 

are the specific acoustic admittances of the walls perpendicular to 
the y and z axes, respectively. 

In a cartesian coordinate system the solution to equation (1) 
for the acoustic pressure can be written as 


p(x,y,z,t) •I I 


m*o n*o 


(4) 


and m and n are 


, . cos h'^ym \ cos ( "^sn ) 

where ^ b j sin 7 ‘ 

Integer indices for each mode. The cosine function is chosen for 

symmetric modes (even m,n) and the sine function for antisymmetric 

modes (odd m,n). A Is the modal amplitude coefficient and k , 

mn * 

Trq /b and Trq /d are characteristic values, all to be determined 
ym zn 

subsequently. It is assumed that the characteristic functions 

e 

(^ 's) form a complete set. Only then can the acoustic pressure 

mn 

field be expressed in the series form given by equation (4). The 
characteristic functions however do not in general form an orthogonal 
set. The set of functions is orthogonal in the two special -ases of 
either zero medium velocity or perfectly rigid duct walls. The 
nonorthogonality of the impedes the use of analytical methods 

to obtain a simple closed-form solution for acoustic pressure. A 
numerical method of solution is therefore suggested. 
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The numerical scheme Involves the calculation of three fundamental 

groups of quantities: a) the characteristic values » b) the modal 

amplitude coefficients and c) the total acoustic pressure. In order 

to determine the characteristic values for each mode, equation (4) is 

substituted into the boundary conditions (3) and into the homogeneous 

counterpart of the governing equation (1). The resulting nonlinear 

system of algebraic equations is solved numerically for k^, 

and vq /d by means of the Newton^Raphson iteration technique, 
zn 

Analytic expressions derived for "slightly" soft walls are used as 

initial guesses in order to compute the exact nodal characteristic 

values for any given values of 6 , 6 , and M. 

y z 

The amplitude coefficients downstream (A and upstream 

(A** ) of the source, are calculated by substituting equation (4) into 

mn 

the governing equation (1). This equation is then Integrated across 
and along the duct. Application of pressure and pressure gradient 
discontinuity conditions at the source location in the resulting equa- 
tion gives two sets of linear algebraic equations with an infinite 

number of unknowns (A and A ). This sytem of equations can 

nil nm 

be truncated at a number of terms (modes) which insures convergence 


ox the solution. 
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is then solved 


numerically for the amplitude coefficients. 

Once the characteristic values and amplitude coefficients for 
each mode upstream and downstream are found, the modal pressure can 
immediately be determined. Finally, summation of the modal pressures 


according to equation (4) produces the total acoustic pressure at any 


location in the duct. 


The input quantities needed for the above calculations are the duet 
dinensions, the wall specific acoustic admittances » the flow Mach 
number » the source location* and the driving frequency. The quantities 
printed out by the computer program may include the characteristic values 
and functions* modal amplitudes, modal pressures and total mean-square 
pressure at any point inside the duct. The computer program at its 
present stage is limited by core size to low and medium frequency calcu- 
lations. A typical run using 35 modes requires less than 30 seconds of 
central processing time. 

Several aspects of the acoustic field in a duct with flow have been 
laves tlgated using this solution capability. Comparable solutions do 
not exist in the literature except for the two simple cases which yield 
orthogonal sets of characteristic functions. As a result* the investi- 
gations have been of an exploratory nature examining the variation of 
the mean-square pressure in the duct with frequency* source position* 
and wall Impedance. The present capability of the computer allows 
approximately 30 modes to be summed. For a duct with a cross section 
of 2 meters by 3 meters this allows a maximum frequency of approximately 
250 Hz to be investigated. Hence* any results would find applicability 
to measuresents in the near-field of an acoustic source. 

Based on results from a limited number of test cases* the following 
qualitative statements can be made: 

1) For axial distances greater than 2-3 wavelengths from the 
source* medium motion reduces the acoustic pressure upstream and enhances 
it downstream when compared to the no flovr case. 

2) For hard walls the opposite Is true, i.e., the acoustic pres- 
sure upstream is higher than that for the no flow case, while acoustic 
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pressure downstreao is lower, as expected. 

3) Acoustic pressure levels throughout a duct with hard walls are 
higher than levels In a duct whose walls have finite adnlttance. This 
holds for the particular wall admittance tested and for axial distances 
greater than one wavelength. 

A) An Increase In the source driving frequency Initially Increases 
the acoustic pressure levels across the duct. 

5) The acoustic field in the duct Is sensitive to source location. 

6) In all cases, the acoustic pressure variations across planes 
perpendicular to the duct axis are more pronounced close to the source 
and damp out away from it. 

The above conclusions are drawn on the basis of test runs for a 2 

meter by 3 meter duct with M ■ 0, .3 and P - 0, .1. The source driven 

at 100 Hz and 200 Hz was located at y • z^ - 0 (center of the duct) 

o o 

and y ■ .30m, z ■ .35m. 
o o 

C. High Frequency Theoretical Solution 

In this section a simple theoretical model Is described for jet 
noise measurements In a wind tunnel. For such measurements sound from 
moving acoustic sources (l.e. , turbulent eddies) Is measured by micro- 
phones located In a moving acoustic medium. A simple theoretical model 
for this situation is an acoustic point mass source which moves with 
constant speed through a uniform acoustic medium moving In the same 
direction with constant speed. This model is used to calculate 
explicitly the far-field acoustic pressure, the Doppler shift in fre- 
quency, and the acoustic power output of the source. Since the effects 
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of tha wind tunnel walla are absent from this model, the results are 
applicable to high frequencies where far-field measurements may be 
obtained within the "hall radius" of the wind tunnel. 

The convective wave equation which describes the acoustic field 
of a mass source in a uniformly moving acoustic medium has been given 
previously by equation (1) in Section B. k point mass source which 
9K>ves with constant speed V In the same direction as the medium is 
specified by expressing the source strength as follows: 

q<x,y,z,t) - q(t)6(x-Vt)fi(y)6(z) . (5) 

Since there are no boundaries in this model, the convective wave 
equation with the source strength specified above is solved for out- 
going waves. The solution for the acoustic pressure is obtained by a 
standard method using a series of coordinate transformations. In order 
to correspond with "in wind" measurements of jet noise, the expression 
for the acoustic pressure is expressed In terms of the distance and 
angle between the "In wind" observer and the source at the time that 
the source emits the signal. The "Doppler shifted" frequency which Is 
observed owing to the relative motion of the source and the observer 
is also calculated from the solution. 

The relationship between far-field measurements of the acoustic 
pressure and the total power output of the source may be obtained by 
using the expression for the acoustic intensity. Two different expres- 
sions appear In the literature for the acoustic Intensity for a 
uniformly moving acoustic medium. In order to ascertain the appropriate- 
ness of either of these expressions to the problem under consideration. 
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both were ueed to obtain the corresponding results for the acoustic 
power. 

The results calculated fron this nodal Indicate a relatively 
strong Influence of the nediun notion on the accustlc field of a moving 
point source. The directional dependence of the mean-square acoustic 
pressure in the far field is easily calculated from the solution. For 
a constant source Mach number* the mean-square pressure decreases for 
all observer locations as the value of the medium Mach number Increases. 
This occurs even though the source strength remains unchanged. If the 
mediu«< Mach number is different from zero* the mean-square pressure 
measured at an angle of ninety degrees to the source direction retains 
a dependence on both source and medium Mach numbers. 

For a constant source Mach number, the ratio of observed to 
emitted frequency Is found to decrease for all observer posltloiiS as 
the medium Mach number is Increased. There is a Doppler shift in fre- 
quency found at an angle of ninety degrees to the source direction 
when the medium is in notion. 

Of the two expressions for the acoustic power which correspond 
with the two definitions of the acoustic intensity* one expression 
contains results which are unreasonable. The second expression gives 
the acoustic power output as a function of both the relative speed 
between the source and medium and the medium speed. This suggests 
that only the latter expression for the acoustic intensily it accept- 
able for the model calculation and perhaps for any free-fleld 


calculation. 
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D» Expcrlnental Investigation with Analytical Conparlson 

Aspects of sound propagation In the 2.13 neter by 3.05 neter wind 
tunnel of NASA Ames Research Center have been Investigated experimentally. 
The primary objective of these experiments has been to explore the effects 
of a mean flow on sound propagation In the wind tunnel. To meet this 
objective, the propagation of acoustic pulses was investigated. 

The acoustic source for these experiments is an AR-1 woofer. This 
speaker Is mounted flush against a plywood window near the middle portion 
of the test section of the wind tunnel. A small hole is drilled through 
the plywood near the center of the speaker. A gated sinusoidal pulse is 
obtained by passing the output of an oscillator through a General Radio 
model 1396-B tone-burst generator. This device Is set to provide a 
pulse of four cycles of the oscillator frequency. The pulse repeats 
Itself typically every one-half second. The oscillator frequency for 
these tests ranged from 50 to 150 Herts. The signal from the tone-burst 
generator Is amplified by a power amplifier which drives the speaker. 

The acoustic signal Is received by two BrOel and Kjaer one-half 
Inch condenser microphones. Both are flush-mounted In the ceiling of 
the wind tunnel. One Is positioned approximately one meter upstream 
from the source and the second is positioned at an equal distance 
downstream. 

When the wind tunnel Is operating, there Is a relatively large 
level of background noise received by the microphones. The primary 
sources of this noise are the wind tunnel drive fan and motor and the 
turbulent boundary-layer along the walls of the wind tunnel. The 
amplitude of the background noise at the higher tunnel speeds is far 
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in access of the maxlimm pulse amplitude which can be generated by the 
equipment used. The pulse signal however can be "recovered" from the 
noise by suitably averaging the signal over many realizations of the 
pulse. The microphone outputs are therefore connected to a SAICOR 
model SAX 43A correlation and probability analyzer which Is used In the 
signal-enhancement mode. After many pulses are averaged, the resulting 
signal is plotted on paper by an X-Y recorder driven by the analyzer. 

All of the pulse shapes recorded In this manner have certain 
similar characteristics. The initial four cycles of the signal are 
generally Increasing In amplitude. The r«&alnlng portion of the signal 
(which will be referred to as the "tall") diminishes In amplitude with 
time. The specific shape of the tall Is strongly dependent on the 
oscillator frequency* The effect of flow on sound propagation Is examined 
by comparing the pulses received upstream and downstream of the source. 

The amplitude of the signal received upwind uf the source is greater 
chan that received downwind, as expected from theoretical considerations. 
The magnitude of the* ratio of amplitudes received upwind to downwind 
depends on the oscillator frequency. There are also differences in the 
arrival time of the pulses upstream and downstream, as expected. 

As a theoretical model of this experlm^t pulse propagation is 
onsidered from a point source located In a rectangular duct of infinite 
length. The duct has a constant cross section and contains a uniform 
flow with velocity V in the positive x-direction. The walls of the 
duct are perfectly rigid. Although not faithful in all respects to the 
experimental conditions, this model does contain many of the basic 
elements encountered during the wind tunnel tests. The governing equa- 
tion which describes sound propagation from a point source in a uniformly 
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movlng acoustic nedlum Is given by equation (1) in Section B. For a 
point mass source located at position 

q(o,y^,z^,t) - q(t)6(x)5(y-y )6(z-z ) . (6) 

o o o o 

The boundary conditions imposed by the rigid walls are 

“ 0, at the walli^. (7) 

3y 3z * 

The problem of pulse propagation in the duct is formulated by 
supplementing equations (1), (6) and (7) with suitable initial con- 
ditions. Before the pulse is initiated at time zero, we assume that 
there is no acoustic disturbance present in the duct. In order to 
compare with the experimental data, we take the source strength to be 
a four-cycle sinusoidal pulse defined as follows 

1 0 t<0 

B sin(2irf t) 0<t<4/f (8) 

o o 

0 r.>4/f 

o 

The governing equation is solved using the Laplace transform. The 
acoustic pressure as a function of time is obtained from this solution 
by numerically evaluating the inverse Laplace transform. This pro- 
cedure permits the theoretical results of pulse propagation in a duct 
with flow to be calculated in a straightforv/ard manner for a realistic 
input pressure excitation. As perusual, pulses with significant high- 
frequency content require longer computation time than pulses of lov; 
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frequency. The results of these calculations are found to be In very 
good qualitative agreement with the experimental data. Quantitative 
differences are found in the amplltudr. of the tail of the pulse. Such 
deviations from experimental results performed in a wind tunnel are to 
be expected from the use of a model based on a duct id.th uniform cross 
section. Nevertheless* pulse propagation in the wind tunnel appears to 
be well described by the theoretical model for the cases examined. 


